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Abstract. The coupled states A*E+ ( 1 D + 1 S), c 3 n u ( 3 P + 1 S) and a 3 E+ ( 3 P + 1 S) of the calcium 
dimer are investigated in a laser induced fluorescence experiment combined with high-resolution Fourier- 
transform spectroscopy. A global deperturbation analysis of the observed levels, considering a model, 
which is complete within the subspace of relevant neighboring states, is performed using the Fourier Grid 
Hamiltonian method. We determine the potential energy curve of the A 1 E+ and c 3 n u states and the 
strengths of the couplings between them. The c 3 n u and a 3 E+ states are of particular importance for the 
description of collisional processes between calcium atoms in the ground state 1 So and excited state 3 Pi 
applied in studies for establishing an optical frequency standard with Ca. 

PACS. 31.50.Df Potential energy surfaces for excited electronic states - 33.20.-t Molecular spectra - 
34.20.Cf Interatomic potentials and forces 



1 Introduction 

It was recently demonstrated that an optical frequency 
standard using an ensemble of ultra- cold calcium atoms, 
which is probed on the intercombination transition 3 Pi 
(m.j = 0) «— 1 So by a high precision laser, has the poten- 
tial to exceed the microwave cesium clock in stability and 
accuracy ^ ne performance of such an optical clock 

depends on the possibility to reduce or to correct for the 
effects of any disturbances during the laser interrogation 
of the transition. Collisions between calcium atoms lead 
to a frequency shift of the transition, which limits the 
achievable accuracy. Therefore, the investigation of the 
collisional processes between atoms in the involved 1 So 
and 3 Pi states is of particular importance for the improve- 
ment of the accuracy of the frequency standard. Such high 
accuracy is required for testing fundamental theories for 
the quest of combining general relativity with quantum 
mechanics or to search for cosmological variation of natu- 
ral constants. 

In this perspective we have reported previously on the 
study of the interactions between two 1 So calcium atoms 
|!J1|4"|. In this article we present our investigations on the 
interactions between one atom in the ground state and one 
excited atom in the 3 Pi,2 state. With the knowledge of the 
molecular potentials correlated to the 3 Pi + 1 So dissocia- 
tion limit, we would be able to calculate photoassociation 
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spectra near the intercombination line and might gain a 
better understanding of trap losses occurring during the 
cooling cycle on this transition. 

Like for the atomic intercombination line, the observa- 
tion of the c 3 IT u and a 3 £+ states by a direct excitation 
from the ground state is rather inefficient. However, the 
c 3 n u state is strongly coupled via spin-orbit interaction 
to the A 1 S+ state, which is reachable from the ground 
state. The a 3 E+ is indirectly coupled to the A : S+ state 
by spin-orbit and rotational coupling via the c 3 n u state. 
Therefore, we can investigate the A state and observe its 
perturbed level structure in order to obtain information 
about all potentials of the coupled states and the cou- 
pling strengths. Generally, we are observing also levels of 
the triplet states taking advantage of their singlet charac- 
ter. 

Therefore, we aim to determine the potential energy curves 
of the coupled-states system A X E+ ( X D + 1 S), c 3 n u ( 3 P 
+ 1 S) and a 3 E+ ( 3 P + 1 S) (see figured). We will per- 
form a global deperturbation analysis of the observed lev- 
els, which enables us to treat all the perturbations of the 
levels, simultaneously. This method will provide physical 
parameters for the potentials and the coupling strengths 
while local deperturbation methods will only give phe- 
nomenological parameters |3j . We collect spectroscopic in- 
formation of a whole set of rovibrational levels of the cou- 
pled states. The shape of the potential curves will be ad- 
justed by a fitting procedure. 
This article is organized as follows. In section [2 we re- 
view the previous investigations of the coupled states. The 
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Fig. 1. Simplified potential scheme of the calcium dimer. The 
potential curves of the A X E^" and c 3 Ii u states are determined 
in this study. The curve of the X 1 S^" state was published in [1] 
and the curve of the B 1 E+ state is given in HJ. The potential 
curve of the a 3 Ej , 3 A U , 3 II U states have been obtained by 
Czuchaj et al. ,T2l by ab initio method. 

experimental set-up and the spectroscopic observations re- 
alized by laser induced fluorescence and filtered laser ex- 
citation spectroscopy are reported in section We derive 
the Hamiltonian of the coupled states in section 21 and 
describe, in sectional the global deperturbation analysis 
for determining the potential energy curves of the coupled 
states and coupling interactions. In section H3 we discuss 
the results of the deperturbation analysis and conclude 
in section by a comparison with earlier results and by 
proposing further studies. 



2 Previous experiments of the coupled system 

The A 1 E+ — X 1 !]^ system was observed by Bondybey 
and English in a supersonic jet created from the vapor- 
ization of calcium metal by a pulsed laser. The spectro- 
scopic interrogation was realized by a pulsed dye laser sys- 
tem limiting the experimental uncertainty to a few tenths 
of a wave number. Several bands starting from low rovi- 
brational levels of the ground state were observed. 
Hofmann and Harris published the results of a more sys- 
tematic investigation of the A 1 E+ - X : E+ system [H]05]. 
They performed laser induced fluorescence spectroscopy 
in a heat pipe oven, and applied the filtered laser excita- 
tion technique (FLE) to observed band heads. They stated 
an experimental uncertainty of 0.01 cm -1 . From the ob- 
servation of P-R doublets they assigned the progression 
to the A 1 ^^ — X 1 !]^ system, and measured 720 lines 
corresponding to 340 different rovibrational levels of the 
coupled states A and c. Their observed and assigned tran- 
sitions are available through the original publication 9 , 



therefore we used their data in our own analysis. Their 
assignment of the lowest observed vibrational level of the 
A state is vh&h = 7. This high lying number 7 leaves 
room for different assignments. Therefore, vibrational as- 
signment will be re-established using the additional and 
more precise spectroscopic data we have measured. This 
notation v H&H for the vibrational quantum number v will 
be used here on to distinguish the earlier from a new as- 
signment given below. 



3 Experimental set-up and spectral 
observations 

The calcium dimer was formed in a stainless steel heat 
pipe oven, characteristics of which are described in a pre- 
vious publication JJj. The pipe was filled with approxi- 
mately 5g of 40 Ca (99, 5% purity) and with argon as buffer 
gas at a pressure around 50 mbar. We operated the oven 
at temperatures between 1240 K and 1275 K. The sam- 
ple was irradiated along the pipe axis by different laser 
sources. A frequency stabilized, linear dye laser (Coher- 
ent 599) pumped by an argon ion laser (Innova 400), us- 
ing Pyridine 1 and DCM dyes, was run in single mode 
with a typical output power of 70 mW. Transitions of the 
A 1 E+ - X X E+ system in the intervals from 13850 cm -1 
to 14550 cm^ 1 and from 14900 cm^ 1 to 15650 cm -1 are 
excited. We collected the induced fluorescence through a 
Fourier-transform spectrometer (Bruker IFS 120HR) by a 
broad band photomultiplier. This total collection results 
to the LIF data set. Each recorded spectrum exhibits a 
single progression of rotational doublets as expected from 
the selection rule A J = ±1 for the rotational quantum 
number J for E <-> E transitions. Few spectra show sev- 
eral progressions due to the overlap of exciting transitions 
with the laser frequency within their Doppler linewidth. 
Because the rotational spacing of the doublets is deter- 
mined by the ground state and thus well known 0], the 
rotational assignment of the observed lines is obtained 
without ambiguity. The strongest lines are surrounded by 
collisionally induced satellites. These satellite lines permit 
to measure the position of neighboring rovibrational levels 
around the excited levels. Figure[2]shows the P-R doublets 
progression from the excitation transition (vh&h = 12, J 
= 29)^(v"=0, J"=28) at 15514.843 cm" 1 . On the insert 
a zoom around the fluorescence lines P(28) and R(30) of 
vx= shows the relatively large number of collision in- 
duced satellites (24 satellites corresponding to 12 different 
J values). 

Due to the possible excitation of different classes of ve- 
locities depending on the frequency of the laser, a Doppler 
shift can occur. Before recording with the Fourier spec- 
trometer, we ensured that the excitation was tuned to 
the maximum of emission for the selected A-X transi- 
tions for which we directed the fluorescence, emitted in 
the forward direction of the laser, to a 1 m monochroma- 
tor (GCA/McPherson Instruments). It was used with a 
band pass of about 2 to 5 cm -1 width. The center of the 
frequency window was set to the position of a line of the 
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Fig. 2. Vibrational progression obtained from the excited 
(vh&h = 12, 29)<-(0, 28) transition of A X E+ - X*E+ . 



induced fluorescence progression with favorable Franck- 
Condon factors, involving typically a low lying vibrational 
level of the ground state. The light passing through the 
monochromator was detected using a broad band Hama- 
matsu photomultiplier (R928). The detection of the fluo- 
rescence intensity allowed setting and fixing the laser fre- 
quency to the maximum of emission. 
In this way, fluorescence progressions have been observed 
from the systematic excitation of rotational levels of vh&h = 
1, 2, 3, 6, 7, 8, 9, and 11. Additional rotational levels of 
vh&h = 5, 10, 12, and 13 have been observed. The levels 
below vh&h = 7 are levels, which were not observed by 
Hofmann and Harris UJ. 

Special care should be put on the experimental uncertain- 
ties. By setting the laser to the center of the line, we re- 
duced the possible Doppler shifts to a magnitude of about 
0.006 cm" 1 w 180 MHz. The drift of the frequency of the 
stabilized laser was less than 10 MHz per hour, which is 
sufficiently smaller than the desired setting of the cen- 
ter frequency given above and certainly also smaller than 
the Doppler width of the lines in this frequency region 
(~ 1.3 GHz w 0.043 cm -1 ), to allow a stable excitation 
during the time of recording (~ 20 min corresponding to 
20 scans). The resolution of the Fourier-transform inter- 
ferometer was chosen to be 0.05 cm -1 . With a triangular 
apodization, the instrument line width is 0.05 cm -1 . Since 
in single mode operation of the laser, only a selected ve- 
locity class is excited, no Doppler broadening is expected, 
and velocity changing collisions play no role under these 
conditions. The lifetime of the A state is 57±5 ns giv- 
ing a homogeneous broadening of about 6 xl0~ 4 cm -1 , 
which is negligible compared to the instrumental broad- 
ening. The broadening due to the size of the aperture of 
the instrument (~ 1.3 mm) is in the order of 0.025 cm -1 . 
Taking the considered effects together, we expect a line 
width of 0.056 cm -1 . The measured widths of the lines 
with good signal-to-noise ratio (SNR), i.e. higher than 10, 
were ~ 0.062 cm . Collisional broadening due to the high 
temperature and the buffer gas is probably responsible for 
the slightly larger value. We have evaluated the relative 



frequency uncertainties of the observed lines to be 10 times 
smaller than their line widths. For lines with lower SNR, 
we estimated a higher uncertainty proportional to the re- 
duced SNR. 

In the frequency range accessible with our dye lasers the 
excitation and detection of fluorescence progressions from 
vh&h = 4 have less favorable Franck-Condon factors (4 
times smaller than for vh&h = 3). In consequence the 
number of collisional satellites with good signal-to-noise 
ratio was reduced to very few. In this case, the Fourier- 
transform spectroscopy does not help to observe a wide 
portion of the vibrational band by the satellites. 
We applied the filtered laser excitation spectroscopy (FLE) 
by recording the emission of Ca2 through the monochro- 
mator. The temperature of the oven was chosen to 1240 K 
since the enhancement of the collisional transfer of popula- 
tion is not needed for the FLE spectroscopy. For absolute 
and relative frequency calibration, the absorption lines of 
iodine from a 60 cm long cell (heated up to 870 K), and 
marker cavity peaks, with 149.7 MHz spectral spacing, 
were recorded. A precise description of the technique is 
given in T . Rotational levels from J = 27 to 77 were ob- 
served for vh&h = 4 in this way. The (4, 61) «— (11, 62) 
measured Ca2 line is presented in the figure|3as example. 
We used the IodineSpec software program JU| to calibrate 
the spectra. It provides a prediction of I2 transitions with 
accuracy better than 25 MHz in this spectral region. The 
width of the calcium lines was 0.042(1) cm -1 correspond- 
ing to the Doppler width of 0.04 cm -1 and a residual 
broadening of about 0.012(3) cm -1 . The width of the io- 
dine lines was in the order of 0.04 cm -1 depending on the 
unresolved hyperfine structure. The absolute experimental 
uncertainty is determined by the precision of the estima- 
tion of the I2 and Ca2 line centers. The signal-to-noise 
ratio was sufficiently high for the I2 and Ca2 lines to allow 
a determination of their centers better than 0.004 cm -1 
for each of them. The final accuracy of the line position 
was then estimated to be 0.006 cm -1 . 
We have constructed term energies from our observed 
transitions and those reported by Hofmann and Harris [H] 
in adding the term energies of the involved levels of the 
well-known ground state potential 0] . From the LIF data 
set, the term energy of a level of the A state is obtained 
several times since the observed spectrum consists of a 
rovibrational progression from that selected upper level. 
Within each observed progression we have averaged the 
constructed term energies weighted by their relative un- 
certainties. To the resulting uncertainty from this average, 
we added the absolute uncertainty 0.009 cm , which is 
the absolute uncertainty of the frequency given by the 
Fourier-transform spectrometer, to obtain the experimen- 
tal uncertainty of one term energy. Then the term ener- 
gies of the same levels from different recordings have been 
averaged, weighted with their respective experimental un- 
certainty. Finally, the total uncertainty varies from 0.009 
cm -1 to 0.0380 cm -1 . The data of Hofmann and Harris 
consist of P and R transitions for 340 levels of the A state. 
We averaged the two values obtained for each upper level. 
We mention here that the levels observed by the filtered 
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Fig. 3. (4, 61) <— (11, 62) Ca2 calcium line recorded with I2 
absorption lines (reversed here) and peaks of a marker cavity 
FPL The traces of Ca2 and I2 have been smoothed. The fluo- 
rescence detection window was set to the transition (4, 61) — » 
(3, 62) at 14284.9 cm _1 and a width of 8v = 2cm -1 . 



laser excitation technique have an experimental uncer- 
tainty smaller than those obtained by LIF and those of 
Hofmann and Harris' data set. Only level vh&h = 4 has 
been observed with this technique. In order to avoid a too 
high weight for the fits on this single vibrational level com- 
pared to all the others, we have increased their uncertain- 
ties to 0.015 cm which is equal to the typical uncertainty 
of the LIF line with good signal-to-noise ratio. 
Figure 0] shows the obtained term energies with respect 
to the minimum of the ground state potential 4 . The 
progressions as function of J(J+1) for each vibrational 
levels of the A state are clearly visible. In the insert, the 
term energies of vh&h = 9, 10, and 11 around J = 45, i.e 
J(J+1) = 2070, are presented as examples, to show the 
perturbed rotational structure of the levels of the A 1 S+ 
state. Around a crossing with a vibrational level of the 
c 3 II u state, the mixing between the two states becomes 
important. The levels of the triplet c state get sufficient 
singlet character to become observable yielding extra lines 
in the spectrum. The energy positions of such levels are 
important since they carry a lot of information about the 
c state and the coupling strength. 

We have compared the data set from Hofmann and Har- 
ris to our own measurements by calculating the differences 
between term energies formed with their and our observed 
transitions. The overlapping set consists of 64 common 
levels. The data set from Hofmann and Harris is shifted 
in average by 0.037 cm -1 to higher energy compared to 
our measured set. The standard deviation of this shift is 
0.038 cm -1 . This spread is larger than our experimen- 
tal accuracy and the one claimed by Hofmann and Harris 
(0.01 cm^ 1 ). We assign to all their data an error of 0.038 
cm -1 . Later on, the global treatment shows that the data 
set of Hofmann and Harris should be lowered by 0.028 
cm^ 1 to give the best fit. We remove from the data set of 
Hofmann and Harris the 64 common levels and keep ours, 
since they have been observed with a higher accuracy. The 
total data set consists of 502 term energies. 



The LIF and FLE spectroscopic techniques enabled us to 
observe six lower vibrational levels of the A state com- 
pared to the data set observed by Hofmann and Harris. 
Local perturbations have been identified at lower energy 
than the energy of the lowest perturbation observed in 
their study, which they proposed to be caused by v c = 
of c 3 n u (0 + ). Therefore, these observations show already 
the necessity of a revision of their vibrational assignment 
for the c 3 n u state. Since more information on the lower 
vibrational levels of the A 1 E+ has been collected, the vi- 
brational assignment of this state should be examined, 
too. A list of all known transition frequencies of the sys- 
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Fig. 4. Term energies with respect to the minimum of the 
ground state potential 4 obtained in this study and in the 
study of Hofmann and Harris. A zoom of the region delimited 
by a square is presented in the insert. Extra levels are clearly 
visible for several J values, for instance for the rotational ladder 
of vh&h = 9, 10 and 11 shown in the insert. 



tern (A X S+ , c 3 n u ) — X 1 !]^ and the derived term energies 
are available in the supplementary Online material. 



3.1 Vibrational and rotational assignments 

From the first inspection of the data we assumed that no 
local perturbation was found for v h&h < 4, which was 
later confirmed by the complete analysis. The coupling 
between the A state and the c state affects the vibrational 
levels below vh&h = 4 only by a global shift of the rovibra- 
tional levels towards lower energies (see [2] ) ■ The variation 
of the shifts with vh&h and J is monotonic compared to 
the local perturbation at higher vn&H-values. These shifts 
are also much smaller than the vibrational spacing of the 
levels. We can thus consider that although a one-channel 
fit of such levels will not provide a good reproduction of 
energy level positions, it will nevertheless be sufficient for a 
good estimate of Franck-Condon factors (FCF). We tested 
different assignments of the A state by comparing the in- 
tensity pattern of the fluorescence progressions to the FCF 
calculated from the fitted potential corresponding to each 
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Fig. 5. Calculated FCF for different assignment of the A state 
levels, compared to the observed fluorescence progressions ob- 
tained by the excitation of the transitions (vh&h = 1, 61) <— 
(8, 62), (vh&h = 2, 57) <- (6, 58) and (v h &h = 3, 47) <- (4, 
46). Agreement is obtained for Av= +1. 



assignment. We tested the following assignments Av = 
— 1, Av = 0, Av = +1 and Av = +2 with Av = v now — 
vh&h- The fluorescence progressions obtained by the exci- 
tation of the transitions (vh&h = 1, 61) ^— (8, 62), (vh&h 
= 2, 57) <- (6, 58) and (v H&H = 3, 47) <- (4, 46) are 
compared to the predicted Franck-Condon factors for the 
different assignments, and normalized to the most intense 
line of each progression. These comparisons are presented 
in figure In this figure the assignment corresponding 
to the Av = +2 has been omitted for saving space. Con- 
vincing agreement is obtained for Av — +1. Therefore all 
observed vibrational levels have been shifted by +1 for 
the A state. This correct assignment will be used in the 
remainder of this article as va = vh&h + 1< 
To determine the vibrational assignment of the c state, 
we have to find the position of the lowest local pertur- 
bation on the rotational progressions of the vibrational 
levels of the A state. This perturbation is then assigned 
to a crossing with v c — of the c 3 ±l u (f2 = 0) state. In fig- 
ure El an estimate of the effective rotational constant B v 
= A(E)/(4J-2) is presented for levels v A = 2 to 8. A(E) is 
the spacing between two observed consecutive rotational 
levels (J — 2) and J. We see on the lower graph that the 
-Revalues of va = 7 and 8 abruptly decrease due to local 
perturbations. Despite that only two -B„-values are ob- 
tained for v A = 6, the significantly lower values clearly in- 
dicate the presence of a local perturbation. If we compare 
the magnitude of the variation of the rotational constants 
for va = 3 and 2 with the higher vibrational levels, effects 
of local perturbations on these two levels should be visi- 
ble in the range of observed rotational numbers, despite 
the observation gaps in the series of levels. But, such a 
local perturbation is not present. Instead we have smooth 
and slight decrease of the effective rotational constant. A 
relatively fast change of the effective rotational constant 
is visible for va = 5 at J values approaching 77. A mod- 
erate decrease of the i3„-values for va = 4 is observed at 
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high J indicating the onset of a local perturbation for this 
vibrational level. Considering that no local perturbations 
are present for va = 3 and 2 for J < 109 and assuming 
that the same vibrational level of the c state causes the 
two local perturbations on va = 4 and 5, we conclude that 
these two local perturbations are the lowest in energy, and 
are caused by the level v^ =0 = of the c 3 II u state, which 
is the lowest in the whole manifold of 3 IT U . These obser- 
vations and assumptions will lead to a consistent picture 
of the perturbations as we will show in section [5] 
Even if there are no crossings with the c state levels of 
the lower levels of the A state they are globally shifted to 
lower energy due to their couplings to all the levels above. 
Consequently, a global treatment of the coupled states is 
required. 

Figurc0gives the range of vibrational and rotational quan- 
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The blue open circles are the data from Hofmann and Harris 
[3- Levels corresponding to extra lines have received the same 
vibrational quantum numbers as those of the A state levels. 

turn numbers observed in this work and those reported by 
Hofmann and Harris. 



c 3 n u and a 3 E+ both correlated to the (4s4s 1 S + 4s4p 3 P) 
asymptote, see figure ^ Other states dissociating at the 
asymptote 4s4s 1 S + 4s3d 3 D can couple to the considered 
states (see fig. but they are expected to lie above the 
triplet manifold according to the ab initio calculations of 
|12j. Their coupling will be much less pronounced due 
to the absence of resonance effects, and will not be in- 
cluded in our analysis. Consequently we will derive poten- 
tials that might contain these small influences, and, in this 
respect, should be considered as effective potentials. Two 
potentials, the 3 n„ ( X S + 3 D ) and the 3 A U (*S + 3 D ) are 
deep enough according to ab initio calculations to overlap 
and cross the repulsive branch of the c state (see fig. . 
But this overlap appears at the very bottom of these po- 
tentials, where the density of levels is small. The magni- 
tude of the perturbations will be most likely small, because 
of low Franck-Condon factors between these states and the 
states A and c. The influence of these states is neglected in 
our deperturbation analysis. For a similar reason the in- 
fluence of the a 3 E+ state is small since it is much deeper 
than the c state to which it is coupled. Near the asymp- 
totic limit ( 3 P + 1 S) where the energy spacing between 
the potential of the c state and the a state becomes smaller 
than the spin-orbit interaction the coupling between these 
states is not negligible. The way to include this state in 
our analysis will be discussed in section !?). 31 



4 The Hamiltonian for the coupled A, c and a 

states 

The Hamiltonian of a diatomic molecule in the body-fixed 
frame xyz can be expressed as: 



H = V BO + T N 



jjROT _|_ jjREL 



(1) 



where V BO is the Born-Oppenheimer potential matrix, 
and T N is the nuclear radial kinetic energy operator. H ROT 
is the rotational Hamiltonian given by: 



H 



2/iR 2 



[R] 2 

[(J 2 - J 2 ) + (L 2 



(S 2 -S 2 )+ 



(L+S_ + L_S+) - (J+L_ + j_L+) - (j + S_ + J_S+) 



spin — electronic 



L — uncoupling 



S — uncoupling 



(2) 



where J = R + S + Lis the total angular momentum 
(the nuclear spin of calcium I =0), R is the nuclear rota- 
tion operator, S is the total electronic spin, and L is the 
total electronic orbital angular momentum. We use the 
ladder operators 0± = O x ± i O y where O stands for J, 
L and S. H REL contains the spin-orbit interaction H so , 
the spin-spin interaction H ss and the spin-rotation inter- 
action H SR . 

Here we limit ourselves to the coupled manifold of states 
A^^ correlated to the (4s4s 1 S + 4s3d 1 D) asymptote, 



4.1 Basis functions 



In Hund's case (a) the angular momentum basis functions 
are \a, J, S, f2, A, E >, where A and E are respectively the 
projections onto the molecular axis of the total electronic 
orbital angular momentum and of the total electronic spin 
momentum, Q = A + E, and a stands for all other state 
labels (e.g. electronic configuration, g/u symmetry). The 
two indistinguishable calcium atoms are bosons with nu- 
clear spin zero, therefore levels of the ground state exist 
only for even J and have (e) symmetry, where the (j) basis 
functions are defined by the inversion symmetry (±)(— 1) J 
for molecules with even number of electrons |13| . Thus, 
only rotational levels with (e) symmetry and odd J of 
the A X E+ , c 3 n u and a 3 E+ system can be observed via 
the excitation from the ground state X 1 E+ . In addition, 
perturbation can only occur between levels of same (j) 
symmetry. Therefore, in the following we only consider 
the subspace of levels with (e) symmetry. The properly 
symmetrized functions with (e) symmetry and odd values 
of J are the following C[TT]1: 

| 2S + 1 A ,J,e> = 
(2)- 1 /2[| a;J;A ,S, E,0 > -(-l)- s |a, J,-A,S,-E,-0 >] 



(3) 
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Then, the normalized wave functions for the five involved 
states and J odd are: 

\A 1 T,t,J,e> = \a, J, 0, 0, 0, > 
|c 3 n u (0 + ),J,e > = (2)" 1/2 [|a, J, 1,1,-1,0 > +\a, J, -1, 1, 1, >] 
|c 3 n u (l),J,e > = (2)- 1/2 [|a,J,l,l,0,l > +\a, J, -1,1,0,-1 >] 
|c 3 n u (2),J,e > = (2)- 1/2 [ja,J,l,l,l,2 > +|a, J,-l, 1,-1,-2 >] 
|a 3 E^ (1), J, e > = (2)" 1/2 [|a, J, 0, 1, 1, 1 > +]a, J, 0, 1,-1, -1 >] 

(4) 

For simplicity the labels (e) and a will not be repeated 
hereafter. 



it contributes as: 

< J, A, S, E, ft|H JS |J, A, S, E— 1, ft+1 >= 



2^R : 



4.2 Matrix elements 

Matrix elements for the different terms of the Hamiltonian 
of eq.£Q) can be found in |1 1 j and an example of applica- 
tion in the case of K2 in |S] . We give here the matrix ele- 
ments concerning our case first for the diagonal elements 
and second for the off-diagonal elements. 
Diagonal matrix elements of H ROT are: 

< J,A,S,E,ft|H ROT |J,A,S,E,ft > = 
h 2 

x [J(J + l)-ft 2 + S(S + l)-E 2 ] (5) 



2^R 2 



Diagonal matrix elements for the different contributions 

in H REL are the following. 

For the spin-orbit interaction H so : 

< J, A, S, E, ft|H so |J, A, S, E, ft >= AEA(R) (6) 

A(R) is the spin-orbit function to be determined. The ro- 
tational interactions couple states within the same multi- 
plicity. 

For the other terms of the relativistic Hamiltonian H REL 
which contribute with a much smaller magnitude than the 
previous one we have for the spin-rotation Hamiltonian 
ED: 

H SR = 7 (R) R • S = 7(R) (J — L — S) • S (7) 

< J,A,S,E,ft| H SR |J,A,S,E,ft >= 

7 (R) [E 2 -S(S + 1)] (8) 

And for the spin-spin Hamiltonian H ss : 



H 



ss 



e(R) (3S 2 - S 2 ) 



leading to: 



<J,A,S,E,ft| H bb |J,A,S,E,ft>= 
e(R) (3£ 2 -S(S + l)) 



(9) 



(10) 



7(R) and e(R) are unknown quantities in our problem. 
For the off-diagonal matrix elements of H ROT , we distin- 
guish the different contributions of the Hamiltonian. The 
S-uncoupling operator is: 



H 



JS 



2^R : 



■(J+S_ + J_S+) 



(11) 



x v/J(J + 1) - ft(ft + 1) x A /S(S + 1)-E(E-1) 

These contributions are symbolized by for ft = and 

by T X JS for ft = 1. 

The spin-electronic term is: 



H 



LS 



1 



r(L+S. 



2^R 2 

it gives the contribution: 

< J,A,S,E,ft|H LS |J, A ■- l.S.E - [.<?. 

h 2 



(13) 



2^R 2 



x v/S(S + l)-E(E-l) x L(R) (14) 



and is symbolized by Tls for ft = 1 and the function L(R) 
is the expectation value of the L± operator: 

L(R) =< J, A + 1, S, E, ft + 1|L+|J, A, S, E, ft > 

=< J,A,S,E,ft|L_|J,A + l,S,E,ft + 1 > (15) 
The L-uncoupling operator is: 
1 



H 



JL 



2^R 2 



(J+L_ + J_L+) 



(16) 



and gives contributions: 

< J,A,S,E,ft|H JL |J,A- l,S,E,ft + 1 >= 



7 x v /J(J + l)-ft(ft + l) x L(R) (17) 



2/j.R 

which are symbolized by Tj £ for ft = and T l ]L for ft = 1. 
Off-diagonal matrix elements of H REL are coming from 
H so and H SR . From the selection rules (see UU), the spin- 
orbit interaction H so couples the A 1 E+ state to the ft = 
component of the c 3 n u state, and leads to the strong- 
perturbations observed in the rotational energy ladder of 
the A state. Also, the spin-orbit interaction couples the 
a 3 E+ (ft = 1) to the c 3 n u (ft = 1). The contribution are: 

< c 3 n u ((£ ), J|H sc Vr+ J > = X (R) (18) 

< c 3 n u (l u ), J |H so | a 3 E u (l u ), J > = -C(R) (19) 

We have noted the unknown molecular matrix elements 
by x(R) an d C(-R-)- Finally the non-diagonal contribution 
from the spin- rotation Hamiltonian is: 



H SR =( 7 /2)(J + S_+J_S + ) 



(20) 



< J,A,S,E,ft|H SR |J,A,S,E- l,ft + l >= 



7 (i?)\/J(J + 1) - n(n + 1) Vs(s + 1) - e(e - 1) 

(21) 

These contributions are symbolized by T a SR for ft — and 
by for ft = 1. 
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/ H BO 



n c = o 
x(R) 



r LS - C(R) 



fis - C(R) 
-1 + ' + ^- 



In summary, the contributions of the various interac- 
tions are abbreviated by: 



r, 



i 



LS 



2fiR 2 



V2L(R) 



1 JL 



2fiR 2 



1 JS 



2^R 2 



XL(R) r} L 



rj s 



2fiR 2 



2fiR 2 



Vx~^2L{R) 



V2X-A 



7 



s/2X~ 



We have noted X = J(J+1). 

The application to the wavefunctions given in eq. 01 re- 
sults in the 5x5 Hamiltonian matrix (1221) for the con- 
sidered subspace of states, for a given J and (e) symme- 
try. The matrix is ordered by the fi-value of the differ- 
ent states ft s , where s stands for the states A, c and a. 
H?° = V s BO (R) +T N (R) contains the Born-Oppenheimer 
potentials and the kinetic energy of the relative motion. 
The A state potential is represented by V®° . The Sl- 
components of the c state split as Vp(R) = V r c BO (R) + 
(SI — 1)^4(R). The potential of fi = 1 component of the a 
state is described by V^ . 



than half of the smallest local de Broglie wavelength A(R) 
of the relative nuclear motion, following the Nyquist the- 
orem |17| for a proper calculation of the energy positions 
of the levels: AR < nfi/^/2^/XW : where AV is the differ- 
ence between the highest considered energy and the lowest 
minimum of the different potentials. By calling upon the 
Nyquist theorem, we implicitly consider the wave function 
as sine waves. This is an approximation since the variation 
of the amplitude of the wave functions in the classically 
forbidden region is not sinusoidal but exponential. There- 
fore we multiply by a parameter lower than one to set 
AR to lower values than estimated by the above equation 
for getting the uncertainty of the calculated term energies 
smaller than our experimental uncertainty. We found that 
the necessary value for a proper representation should be 
smaller than 0.7 times the Nyquist estimate. We used the 
factor 0.5 to have a correct representation and to avoid an 
unnecessary high number of grid points. The term energies 
are located between 14600 cm -1 and 16272 cm _1 (the zero 
of energy is the minimum of the ground state potential see 
4 ) . The classical vibrational motions are restricted in the 
interval from 3.08 A to 4.52 A for the A state and from 
3.37 A to 8.65 A for the c state in the non-coupled pic- 
ture. For our grid choice ranging from 2.12 A to 11.11 
A the representation of the states of the model (p = 5) 
requires a 1235 x 1235 matrix. 



5 Global deperturbation 

5.1 Fourier Grid representation of the Hamiltonian 

The potentials of the considered sub-space will be deter- 
mined by minimizing the standard deviation between ob- 
served term energies and calculated term energies from 
solving the Schrodinger equation with the Hamiltonian of 
(eq. H22|) ~). The minimization is realized with a non-linear 
fitting routine |14) . The Schrodinger equation is solved 
using the Fourier Grid Hamiltonian method (FGH) [Bill 51 
ITB| . In this method the internuclear distance is discretized 
by a grid with N equidistant points for a length L. In 
this representation the potential and coupling operators, 
which are local in the R-coordinate space, are diagonal. 
In contrast the kinetic energy matrix, which contains the 
derivative with respect to R is non-diagonal. Therefore the 
Hamiltonian matrix of a system of p coupled states has a 
dimension (p-N)x(p-N), composed of p 2 blocks with size 
NxN. 

The spacing between grid points should be at least smaller 



5.2 Construction of potential and coupling functions 

We split the representation of the potentials in three re- 
gions: the repulsive wall (R<R Jlm ), the asymptotic region 
(R>R out ), and the intermediate region in between. For the 
intermediate region, we use the following representation: 

n 

V fit (R)=^a^(R) i (23) 

i=0 

with the analytic function 

« R » - (to) < 24 > 

where {a^}, 6, and R m are free parameters (R m is close 
to the value of the equilibrium distance). We used 13 pa- 
rameters a, ; for the A state and 17 for the c state. The ex- 
perience in our group, obtained by fitting potentials with 
this choice of representation, shows that this manifold of 
parameters should allow correct descriptions of the poten- 
tial curves, for the range of observed vibrational levels and 
for the expected well behaved potentials according to the 
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ab initio calculations from 12 (no shelf or double well 
structure 1 ). We have no direct spectroscopic observation 
for the a 3 E+ state, as it is weakly coupled to the A 1 E+ 
and c 3 n u (SI = 0) states. So, we determine an initial set of 
parameters from the published ab initio curve of ref. |12) . 
These latter parameters were kept fixed during the fit. 
The potentials are continuously extrapolated below Ri nn 
with: 

Vbo(R) = A + B/R 10 (25) 

by adjusting the A and B parameters. Beyond R ou t the 
following forms are considered respectively for the A, c 
and a states: 




V^ (R) = Di 



C 5 A /R 5 



C^/R 6 



(26) 



Intemuclear separation 



Fig. 8. Adiabatic potentials of u-symmetry states correlated 
to the 3 P + X S asymptote. 



V^ (R) = Dl - C 6 C /R 6 - C 8 C /R 8 (27) 

V^ (R) = Dl - C a jR & - C 8 7R 8 (28) 

Rout is different for the three states. The short and long- 
range extrapolation functions are connected near the short- 
est and the largest classical turning points of the observed 
levels. The dissociation limits of the states Df, D c e and 
are calculated with respect to the minimum of the ground 
state potential 0], using its value of the dissociation en- 
ergy Df = 1102.074(9) cm" 1 : 

Df = AE^D- 1 S)+D* 

= 21849.634+ 1102.074= 22951.708(9) cm" 1 and 
D c e = AE^P- 1 S) + D* 

= 15263.003 + 1102.074 = 16365.077(9) cm" 1 

The atomic transition energies are taken from reference 

The parameter C5 was fitted because the A-c coupling in 
the classical forbidden redion of the A 1 E+ state with the 
classically allowed redion of the c gives significant shifts 
of levels. However, our value of C§ should not be consid- 
ered as true long-range coefficent because even the highest 
observed values are far below the atomic asymptote. The 
continuity through the connection point {R^ ut = 4.511 A) 
is ensured by Cg. The Cf and Cg coefficients for the c 3 n u 
and the a 3 E+ states were fixed to the most recent ab ini- 
tio values from ^JJ- The coefficients Cg and Cg were set 
such that the potentials are continuous at the connecting 
points R c out and R^ ut . 



5.3 Four-states model 

Our data set contains only levels, which are accessible via 
the ground state, i.e. which have a strong singlet character 

1 In the ^JCD + 1 S) is lying 2000 cm" 1 above the 3 P 
+ S asymptote. The potential curve of the 1 A U ( 1 D + 1 S) state 
has a similar depth and shape as the A state potential deter- 
mined experimentally. There is probably some wrong ordering 
of electronic symmetry in [12] . 



(> 15% an estimation from the deperturbation analysis 
below). Since the S7 = 1 component of the c 3 II u and a 3 E+ 
states and the 57 = 2 component of the c 3 n u state are 
weakly coupled to the A : E+ state, the number of levels 
of these states which have a significant singlet character 
is low. Hence, the information about the = 1 and 2 
components of the c 3 n u state is low or negligible with our 
data set. Therefore, the information to characterize the 
spin-orbit splitting ^4(R) is limited. We then define the 
fitted potential as the potential of the SI = component 
of the c 3 n u state: 

K fit (R) = V C °(R) (29) 

Consequently the potentials for the other Sl-components 
of the c state are defined as: 

V% od (R) = V c &t (R) + nA(R) (30) 

The a 3 E+ state potential is about four times deeper than 
the investigated energy range below the asymptote 3 P+ 1 S, 
which induces a small grid step in the FGH representation 
(see section l5~T|) . and increases the Hamiltonian matrix 
size in the FGH representation. Since this state is weakly 
coupled to the c 3 n u (SI = 0) and only strongly coupled 
to the c 3 n u (SI = 1) its influence on the description of 
the spectroscopic data is expected to be weak. We took 
into account the coupling between the a 3 E+ (SI = 1) state 
and the c 3 n u (SI = 1) state by replacing the c 3 n(Sl = 1) 
potential by the adiabatic l u (c 3 n) potential with respect 
to the spin-orbit coupling £(R) (fig-Ellj obtained by the 
diagonalization of the simplified 2x2 sub-matrix from eq. 
J221 for each R: 

Sic = 1 o = 1 

fv^ + A(R) r L5 -C(R)\ rrn 
U^-C(R) v a BO ) ^ 

where rotational couplings have been neglected. The model 
is thus reduced to an effective 4-states model, and conse- 
quently the total representation is reduced to a 293 x 293 
matrix. The spin-orbit coupling £(R) has been approxi- 
mated by C(R) = A{R). 

To our knowledge, there are no ab initio calculations avail- 
able, which could properly model the spin-orbit splitting 
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1000 1500 2000 2500 3000 3500 
I(J+1) 



Fig. 9. Observed (full circle) and calculated initial term ener- 
gies for the c state (dotted lines) and for the A state (full lines) 
as function of J(J+1) for the same region as the insert figure 
[I] Vibrational quantum numbers are given at the left for the 
A state and at the right for the c state. 

A(R) and the spin-orbit couplings x(R)- Only the values 
at the atomic limit are known. The spin-orbit splitting 
converges to the atomic spin-orbit value A = 52.940 cm -1 
|18|. The spin-orbit coupling x(R) vanishes for large R 
since the atomic coupling between the (483d) 1 D state and 
the (4s4p) 3 P state is zero because of different parities. We 
have represented these interactions in different R ranges as 
for the potentials (see above). Polynomial functions in the 
inner range allow for sufficient flexibility. For the outer re- 
gion we ensure continuity at the connecting point R c and 
correct atomic values for the asymptote. The representa- 
tion of the functions is expressed as: 

( k(R s ) + k s x (R - R s ) 2 for R < R s 

K (R) = < Ep=o " for R s < R < R c 

[ kI/R 1 + k° for R c < R 

(32) 

where k(R) stands for A(R) or x(R), and asymptotically 
k° — Aq for A(R), and k° — for x(R)- The expansion 
parameter R2 is chosen to be close to the internuclear dis- 
tance corresponding to the crossing between the A 1 E+ 
state and the c 3 II u state, and k p , k s are fitting parame- 
ters. The parabolic extrapolation below R s is used only 
for x(R). ^4(R) is not truncated below R s . 



5.4 Fitting strategy 

The fit of our spectroscopic data to analytical potential 
curves and coupling functions is a highly non-linear prob- 
lem and the parameters to be derived are strongly cor- 
related. Therefore, we should start the fitting procedure 
with the simplest model possible in order to reduce the 
number of parameters and thus the complexity of the 
problem. For this reason, we first limited the model to the 
SI = + components of the A 1 S+ state and the c 3 IT u state 
and simplified the polynomial part of the non-diagonal 



spin-orbit functions (cq. (|32|l ) to an adjustable constant. 
Additionally, the time for achieving convergence of the 
fit depends critically on the initial guesses for the poten- 
tial energy curves and the couplings. Therefore, we con- 
structed initial potentials for the £1 = + components of 
the A 1 !!* and the c 3 IT u states, which yield the best possi- 
ble reproduction of the observed rovibrational ladders. As 
we can see in the figures 0] and El the vibrational and rota- 
tional structure can be well identified for the A state. To 
construct the potentials, we calculated traditional Dun- 
ham parameters from the rovibrational ladder and derived 
RKR potentials with these parameters and converted to 
analytic potentials according to eq. Q23J) by a simple linear 
fit. For the c state, we only have the local perturbations 
to construct an initial potential. The rotational ladder of 
one vibrational level of the c state crosses several rota- 
tional ladders of different vibrational levels of the A state 
at different J- values. This was used to estimate the rota- 
tional constants, and vibrational spacing to obtain an ini- 
tial c state potential. These constructed potentials give the 
rovibrational levels exemplified in figurcEl In this manner, 
we ensured that the ladders of both states cross at the ob- 
served local perturbations. 

We can expect that the coupling of the states leads to a 
shift of each level locally caused by the closest levels of 
the perturbing state and to a global shift coming from the 
accumulated influence of all the other levels. The selection 
of levels following a regular rovibrational series estimated 
from the observed levels, which are not strongly deviating 
from such regular behavior, gives a rough estimation of 
the positions of the levels of the uncoupled states, since 
the global effect is not properly taken into account. In 
a first iteration step, we fitted only the spin-orbit con- 
stant. In this way, we obtained a first guess of the value 
of this constant and the magnitude of the shifts caused by 
the coupling. These shifts have been subtracted from the 
observed term energies, which are not strongly perturbed. 
Then, the initial potential of the A state has been fitted to 
these roughly deperturbed levels. The obtained potential 
has been used again to improve the value of the spin-orbit 
constant. The procedure has been repeated several times 
to improve iteratively the A state potential and the cou- 
pling. Then, this method has been employed once for im- 
proving the c state using this time the strongly perturbed 
levels and the obtained value of the coupling constant. The 
initial potentials and the constant coupling result in a di- 
mcnsionless standard deviation a = 170 between observed 
and calculated levels. The improved potentials obtained 
by the iterative procedure lead to a standard deviation of 
a = 60. 

The obtained potentials and the value of the spin-orbit 
parameter provided the starting conditions for the global 
fit itself. 



5.5 Result with the 4-states model 

Practically, we have extended the fit iteratively, first both 
+ components, second adding f2 = 1 of c 3 II and third 
Q = 2 of c 3 II to avoid to fit directly a large number of 
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free parameters. This is justified since the magnitude of 
the coupling is weak between components with ^ and 
= compared to the spin-orbit coupling between the 
A*E+ state and the c 3 n u (O = 0) state. Therefore, the 
adiabatized = 1 component and the = 2 component u 
of the c state and the diagonal spin-orbit splitting were 
successively introduced when the convergence of the sim- 
plified models was reached. The only contribution, which 
remains neglected in the model, is the spin-spin splitting 
(eq. {HI) since its magnitude is expected to be small com- 
pared to our experimental uncertainty. The parameter e 
was kept at zero. 

Figure shows the quality of the final representation of 
the observed levels, which have been reached with the 4- 
states model. The rotational ladders of va = 7, 12, and 16 
are taken as examples. For each case in the upper part, we 
show the observed levels from which we have subtracted 
the linear function E„(J = 0) + B v x[J(J+l)] in order 
to show better the good reproduction of the perturba- 
tions. For each va , the lower graph shows the differences 
between observed and calculated term energies. We see 
that most of the residuals lie within the experimental er- 
ror bars. The perturbations, which can reach magnitudes 
of 8 cm -1 , are reproduced to tenths of a wave number. 
Therefore, the global reproduction of the position of ob- 
served levels given by the normalized standard deviation 
a = 1.65 is satisfactory corresponding to a standard devi- 
ation of 0.053 cm -1 . This quality is only achievable with a 
model which includes the four considered states. FigurelTD 
shows the improvement beginning with the 2-states model, 
containing only the components O = + of the A 1 £+ and 
c 3 n u states, by the successive addition of the adiabatized 
= 1 component (3-states model) and = 2 component 
(4-states model) of the c 3 II u state. The improvement is 
particularly visible for the A state levels near J = 39 and 
81, which have a clear mixing with the = 1 state and 
a weaker mixing with = 2. Note the different scales of 
the panels. The inclusion of the third and, later on, fourth 
state leads to local improvements but also gives a bet- 
ter overall representation. The relative magnitudes of the 
residuals for the 2-states model compared to the 3-states 
model exemplify clearly the important role of the different 
states. The achieved standard deviation with the 2-states 
model is only a — 4 and a — 1.98 for the 3-states model, 
and finally a — 1.65 for the 4-states model. 
Nevertheless, a completely satisfying reproduction of the 
observed levels within their experimental uncertainties was 
not reached with the considered model. In figures 1101 and 
1111 we see that systematic trends in the residuals remain 
(for instance in fig.^|for va = 16 between J(J+1)« 250 
and 1000 ). We also see that the large residuals are located 
systematically at the position of the strong perturbations 
and they are two to four times larger than the error bars. 
We verified that such trends and large residuals will not 
be reduced by an increase of the number of parameters for 
the potentials. 

Figure El shows the spin-orbit functions obtained by the 
fit. In the range of data (i.e. within the classical turning 
points of observed levels), their variations are slow, despite 
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Fig. 10. The calculated term energies with the 4-states model 
and observed term energies as function of J(J+1) for three se- 
lected vibrational levels va = 7, 12, 16 are presented. A linear 
rotational function E„(J = 0) + B„x[J(J+l)] has been sub- 
tracted from the term energies in order to show the deviations 
from the regular rotational series due to the perturbations. 
The lower graphs present the difference between observed and 
calculated term energies and the experimental error bars. 



the large flexibility that the choice of representation (eq. 
J32| l) offers. Therefore, an increase of the number of pa- 
rameters for the coupling functions will not give a better 
representation of the observations. This can be explained 
by the fact that the data are mainly sensitive to the av- 
eraged strength of the couplings and not so much to the 
details of the shape of the couplings. Hence, obtaining 
more accurate spin-orbit functions from ab initio calcula- 
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Fig. 11. Residuals as function of J(J+1) obtained for the rota- 
tional series of va = 7 with the 2, 3 and 4 - states models. The 
corresponding J quantum numbers at the local perturbations 
are given on top of this figure. 
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Fig. 12. Fitted spin-orbit splitting A(R) of c 3 IT and spin-orbit 
coupling x(R) of A J E+ and c 3 II u (0 + ). 



tions may not help to improve the fit using the present 
data set. 

The derived parameters for the potentials and the cou- 
pling functions are given in tables ^ to 0] The potential of 
the A^^S^" state is determined in the internuclear interval 
from 3.08 A to 4.52 A and the fitted potential for the c 3 II u 
from 3.37 A to 8.60 A. The outer turning points for the 
two potentials are located at relatively short distances cor- 
responding to potential points located 6713 cm -1 below 
the : D + *S asymptote for the A 1 E+ state and 40 cm -1 
below the 3 P + 1 S asymptote for the c 3 II u (0 + ) state. Con- 
sequently, the dispersion coefficients given in tabled to |3] 
should be considered only as extrapolating parameters in 
order to have proper boundary conditions to diagonalize 
the Hamiltonian matrix. 



6 Discussion 

The variations of the residuals seen fig.HOIand fig.llllcould 
be attributed to perturbations caused by other states. But 
they have small magnitudes and appear at positions of 



Table 1. Parameters of the analytic representation of the 
A x Eu state potential. The energy reference is the minimum 
of the ground state potential. Parameters with * are set for 
extrapolation of the potential. 



R < R inn = 3.080 A 


A* 


0.131237254x10 s cm" 1 


B* 


0.23850904xl0 9 cm" 1 A 10 


3.080 A < R < R out = 4.511 A 


b 


-0.57 


Ria 


3.59402020 A 


a 


14106.8528 cm" 1 


ai 


-0.2251508202149158 xlO 2 cm -1 


a.2 


0.1220394308918707x10 s cm" 1 


13 


0.6764477900769186 xlO 4 cm" 1 


fl4 


-0.4630366850561168 xlO 4 cm' 1 


15 


0.1699400434605241x10 s cm" 1 


a 6 


0.9375762656463659x10 s cm" 1 


a 7 


-0.4561172682499647x10 s cm" 1 


ag 


-0.6354314702263585 xlO 6 cm" 1 


ag 


-0.6783215663449778 xlO 6 cm^ 1 


aio 


0.1208346212864403 xlO 7 cm" 1 


an 


0.3797632603225922 xlO 7 cm" 1 


0.12 


0.3059003584697336 xlO 7 cm" 1 


Rout < R 


D£ 


22951.708(9) cm" 1 


C 5 


0.5028185 x10 s cm" 1 A 5 


C 6 * 


-0.1699751X10 9 cm" 1 A 6 


Additional constants: 


equilibrium dist 


ance: R* = 3.595(1) A 


electronic term 


energy: T/= 14106.8(10) cm" 1 



level crossings between the A 1 E+ state and the two com- 
ponents (SI = and 1) of the c 3 II u states. Therefore, possi- 
ble additional perturbing states should only have an indi- 
rect coupling to those latter states or a direct coupling to 
the c 3 II u (O = 1) state. According to the ab initio calcula- 
tions from Czuchaj et al. the bottom of the 3 H U state 
and the 3 A U state, dissociating to 3 D + 1 S, approach re- 
gion of observed levels. That 3 H U state is coupled directly 
to the A state and to all fi-componcnts of the c state. 
Thus, perturbations should also occur at different J and 
not only at J corresponding to rotational level crossings 
between the A state and the c state. Such additional per- 
turbations are not observed within our experimental un- 
certainty. This state is probably less deep than predicted 
and may lead only to a global shift to all the levels. The 
derived potentials in our model could include such global 
shift and would be in this respect effective potentials. 
The 3 A u (f2 = 1) state and the a 3 S+ state are indirectly 
coupled to the = components of the A and c states 
via a direct coupling to the c 3 n u (S! = 1) state by spin- 
orbit interaction. We took the potential curves of these 
states |12j . one by one, and the non-diagonal spin-orbit 
function x(R) to calculate their respective influence on 
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Table 2. Parameters of the analytic representation of the c 3 IT u 
(0 + ) state potential. The energy reference is the minimum of 
the ground state potential. Parameters with * are set for ex- 
trapolation of the potential. 



R < R inn = 3.319 A 


A* 


0.146530549X10 5 cm^ 1 


B* 


0.29313693X10 9 cm" 1 A 10 


3.319 A < R < R out = 8.450 A 


b 


-0.22 


Km 


4.06654500 A 


a 


14841.9738 cm" 1 


ai 


-0.3102692708740178 xlO 1 cm -1 


0.2 


0.1682036434942518 x10 s cm" 1 


0-3 


-0.1609245688758938x10 s cm" 1 


a.4 


0.3095300726978432x10 s cm^ 1 


05 


-0.4053968648443735x10 s cm" 1 


oe> 


-0.2549434419469327X10 7 cm^ 1 


07 


0.2063831175657142x10 s cm" 1 


OS 


0.1659840951559789 x10 s cm" 1 


09 


-0.3493050617277784 xlO 9 cm" 1 


010 


0.3805220943990989 xlO 9 cm" 1 


on 


0.2346473224525490 xlO 10 cm" 1 


012 


-0.5935661971035855 xlO 10 cm" 1 


ai3 


-0.1908247395049612X10 10 cm" 1 


on 


0.2221570121722804X10 11 cm" 1 


015 


-0.2838711585484613 xlO 11 cm" 1 


oie 


0.1186121797412387X10 11 cm" 1 


Rout < R 


Dg(0+) 


16312.139(10) cm" 1 


C 6 from [SI 


0.1186500x10 s cm" 1 A 6 


c 8 * 


0.3427436X10 9 cm" 1 A 8 


Additional constants: 


equilibrium distance: R^=4. 067(1) A 


electronic term 


energy: T e c = 14842.0(10) cm -1 



the position of levels of the coupled states, leading to a 
5-states model. We found that among these states only 
the 3 A„ can create additional deviations to the calculated 
levels. They appear at positions of strong perturbations 
between the levels of the A 1 !]^ and c 3 II u (SI = 0) and 
at positions of weaker perturbations between the A-'-E - }' 
states and the c 3 II u (SI = 1) state. Furthermore, the cal- 
culated magnitude of these perturbations is of the same 
order as the observed larger residuals. The 3 A„ is a good 
candidate to explain why the present model was incom- 
plete for providing a totally satisfactory reproduction of 
the observed levels. For reaching the experimental pre- 
cision by the coupled-states calculation it might be thus 
necessary to extend the 4-states model by the 3 A U (S1 = 1) 
and (SI = 2) states. This implies that spectroscopic data 
for these states are needed, which could be obtained by 
multi-photon spectroscopy or by searching for very weak 
extra lines from the coupling. 



Table 3. Parameters of the analytic representation of the a 3 E u " 
state potential from the ab initio potential published in |12| . 
The energy reference is the minimum of the ground state po- 
tential. Parameters with * are set for extrapolation of the po- 
tential. 



R < Rinn = 2.900 A 


A* 


0.100173775x10 s cm" 1 


B* 


0.25535671xl0 9 cm" 1 A 10 


2.900A < R < R ou t = 8.000 A 


b 


-0.21 


Rm 


3.74686320 A 


a 


8476.274656 cm" 1 


ai 


-0.4989163630417485 xlO 2 cm" 1 


02 


0.5159297957888653x10 s cm" 1 


03 


-0.8701696364839145 xlO 4 cm" 1 


0.4. 


-0.5354588124673624 x10 s cm" 1 


05 


0.9811459564252647X10 4 cm" 1 


a a 


-0.9842985554356317x10 s cm" 1 


07 


-0.6320055597858555 xlO 6 cm" 1 


08 


0.7170441797704500x10 s cm" 1 


ag 


0.2295720046533766 xlO 7 cm" 1 


aio 


-0.2613404819600489xl0 7 cm" 1 


Rout < R 




16365.078(9) cm" 1 


C* 6 from [HO 


0.1855883x10 s cm -1 A 6 


C 8 * 


0.4785476 xlO 9 cm" 1 A 8 


Additional constants: 


equilibrium distance: i?°=3.748 A 


electronic term 


energy: T e a = 8476.3 cm -1 



Notable differences exist with the recent ab initio calcula- 
tions of potential energy curves for these states, in partic- 
ular for the c 3 II u (0 + ) state. The results of the relativistic 
ab initio calculations presented in [20) give a dissociation 
energy of 168 cm -1 , which is in strong disagreement with 
the value of 1470.2(10) cm -1 determined in our experi- 
mental study. Furthermore, the last two potential points 
given in 20; at large internuclear distance lie above the ex- 
pected asymptotic limit, which is in contradiction with the 
attractive long-range behavior of the c 3 II u (0 + ) state (see 
). The dissociation energy of 1575 cm -1 of the c 3 II u 
state derived with non-relativistic calculations by Czuchaj 
et al. 2D is in reasonable agreement (within 6.6%) with 
our experimental value. The agreement with the dissocia- 
tion energy of the A 1 E+ of 9518 cm -1 given in JT2J (Re- 
member the proposed reassignment given in footnote 1) 
compared to our derived experimental value of 8845.0(10) 
cm -1 is in the same order (within 7.6%). 



7 Conclusion 

We have extended the spectroscopy [711311] of the coupled 
system A*E+ (*D + 1 S), c 3 n u ( 3 P + 1 S) and a 3 S+ ( 3 P 
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Table 4. Parameters of the analytic representation of the cou- 
pling functions A(K), C(R) an d x(R) as defined ea. 1321 ea. 1101 
and|Hl 



coupling function 


parameter 




R < R s = 4.149 A 




Xs = -0.756636 x 10 1 cm" 1 A" 2 




R s < R < R c = 8.000 A 




R 2 = 6.00 A 




Xo = 0.364981 x 10 2 cm" 1 


x(R) 


XI = -0.177038 x 10" 1 cm" 1 A" 2 
X2 = 0.186890 x 10" 1 cm" 1 A" 2 




R c < R 




X° = 0.104382 X 10 4 cm" 1 A 




X° = cm -1 




R < R c 




R 2 = 6.00 A 




A = 0.483295 x 10 2 cm" 1 


A(R) = C(R) 


Ai = 0.294419 cm' 1 A" 1 

A 2 = -0.2914325 x 10" 1 cm" 1 A" 2 




R c < R 




A° = -0.118252 x 10 3 cm" 1 A 




A° = 0.529400 X 10 2 cm" 1 from DU 


7 


0.1632 x 10" 2 cm" 1 


e 


fixed to zero 



+ 1 S) using the laser induced fluorescence and the filtered 
laser excitation technique. Rotational series of six new vi- 
brational levels at the bottom of the A state have been 
observed compared to the previous study of Hofmann and 
Harris ■ It allowed us to reassign the A state with great 
confidence. The lowest perturbations of the A state lev- 
els have been identified and consequently the c state has 
been reassigned. We have developed a description of the 
observed levels in which the spin-orbit couplings and the 
rotational couplings between the A 1 S+ , c 3 II u and a 3 E+ 
are modeled. The Fourier Grid Hamiltonian method in- 
cluded in a fitting procedure allowed a global treatment 
of the coupled states, and led to a detailed description 
of the spectroscopic observations compared to the previ- 
ous derivation from Hofmann and Harris UJ. They have 
obtained a reproduction of the observed 340 term ener- 
gies of 0.13 cm -1 by a local deperturbation analysis. 183 
coupling parameters and 16 Dunham coefficients were de- 
termined to reach this quality. These large numbers are 
due to the fact that the coupling of each level to all the 
other levels of the coupled states are modeled with local 
perturbations introducing for each new perturbing level a 
new parameter. In addition, there is some ambiguity in 
the choice of levels to be considered as influencing. This 
shows the limits of local methods for the treatment of rel- 
atively strongly coupled states. 

We have determined potential energy curves for the A 1 E+ 
state and the c 3 n u state as well as spin-orbit coupling 
functions. For the two effective states (A 1 ^ , c 3 n u (0 + )) 
and the coupling functions 39 parameters have been de- 



rived from the fit and 22 parameters fixed for the a 3 E+ 
(1) state potential or for the continuation to small or large 
R . Small deviations remain with the present model. We 
believe that such effects can be explained by the neigh- 
boring 3 A U ( 3 D + 1 S) state, which was not reasonable to 
include in the model. More precise data of the studied sys- 
tem, and in particular for the c 3 n u state and the 3 A U are 
needed. 

In their spectroscopic investigations of Ca2 in a supersonic 
jet Bondybey and English j2j have observed band system 
centered around 15000 cm" 1 belonging to the A X S+ — 
X 1 !]^" system and additional bands that were not identi- 
fied (see fig. 1 in jTj). Our analysis allowed us to correct 
the vibrational assignment of these molecular bands and 
to assigned the previously unidentified bands to the c 3 n u 
— X 1 !!^ system (see [B]). Since levels of the c 3 n u state 
have been directly observed in their jet experiment, the 
spectroscopic study of the c state in a molecular beam 
seems promising. Particularly the observation of high ly- 
ing levels of the c state is of interest for a better character- 
ization of the collision processes between one 3 P and one 
1 S calcium atom. The short-range coupling of the other 
state, except the a 3 Sj state, are predicted to be weak for 
the long-range levels of the c state. In this respect the de- 
scriptions will not involve a more complex model than the 
present one. Therefore, this model should provide a good 
understanding of ultra-cold collisions in a Ca trap where 
the intercombination line of calcium is studied to establish 
this resonance as an optical frequency standard. 
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